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Abstract 

We present a first attempt to experimentally extract an effective strong coupling constant that 
we define to be a low Q 2 extension of a previous definition by S. Brodsky et al. following an initial 
work of G. Grunberg. Using Jefferson Lab data and sum rules, we establish its Q 2 - behavior over the 
complete Q 2 -range. The result is compared to effective coupling constants inferred from different 
processes and to calculations based on Schwinger-Dyson equations, hadron spectroscopy or lattice 
QCD. Although the connection between the experimentally extracted effective coupling constants 
and the calculations is not established it is interesting to note that their behaviors are similar. 

PACS numbers: 12.38. Qk,11.55.Hx 



At experimentally accessible distances, the strong force remains the only interaction that 
resists satisfactory understanding. Quantum Chromodynamics (QCD), the gauge theory of 
the strong force, is well known at short distances (< 10 -16 m) where it is solvable pertur- 
batively. QCD, however, is not perturbatively calculable at larger distances, typically the 
scale of the nucleon radius. Recent precision data on moments of nucleon structure func- 
tions reveal a smooth transition from small to large scales, while in contrast, a 
feature of perturbative QCD (pQCD) is that at Aq C d, the running strong coupling constant 
a s becomes infinite. An approach using an effective coupling constant could reconcile these 
two seemingly paradoxical aspects of strong interaction. 

In lepton scattering a scale at which the target structure is probed is given by the inverse 
of Q 2 , the square of the four- momentum transfered to the target. One way to extract a s 
at large Q 2 is to fit the Q 2 -dependence of the moments of structure functions. Among 
all moments, the Bjorken sum rule is a convenient relation for such an extraction jf|. 
Furthermore, as will be discussed, the Bjorken sum may offer unique advantages to define a 
effective coupling at low Q 2 . 

In the limit where the energy transfer v and Q 2 are infinite, while x = Q 2 / (2 Mis) remains 
finite (M is the nucleon mass), the Bjorken sum rule reads: 

r? _ n _ r; _ r? _ j 1 _ g „ {x)) = (1) 

g\ and g™ are spin structure functions for the proton and neutron. The axial charge of the 
nucleon, g^, is known from neutron /3-decay. for finite Q 2 much greater than Aq CD , the 
Q 2 -dependence of the Bjorken sum rule is given by a double series in Q 2 ~ l (t=2,4... being 
the twist) and in (a s /ir) n . The (a s /7r) n series is given by pQCD evolution equations. At 
leading twist (t = 2) and 3 rd order in a s , in the MS scheme and dimensional regularization, 
we have Q: 

rr* = \qa[1 - - - 3.58 (^V - 20.21 (^)\ (2) 
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The validity of the sum rule is verified at Q 2 =5 GeV 2 to better than 10% [8J]. 

The extraction of the Bjorken integral using data from the Thomas Jefferson National 
Accelerator Facility (JLab) in the Q 2 -range of 0.17-1.10 GeV 2 has been reported recently 
The use of Eq. [2] as an Ansatz for definition of an effective running coupling constant at low 
Q 2 has ambiguities and difficulties as well as some practical advantages. The advantages 



are the following: firstly, it was pointed out in Ref. [10j that the extraction of a s does not 
depend strongly on the low-x extrapolation. Secondly, this flavor non-singlet contribution 
does not mix quark and gluon operators when evolved. Hence the pQCD evolution is known 
to higher order. Thirdly, the JLab data are at constant Q 2 . This avoids a possible ambiguity 
encountered in previous experiments, namely that in order to combine neutron and proton 
data, structure functions must be evolved to a common Q%, which needs a s itself as input. 
This is especially important in our case since we cannot anticipate the value of a s at such 
low Q 2 , due to the breakdown of pQCD. Hence, no evolution of the data to Qq is possible. 
Difficulties within the pQCD approach are: firstly, at low Q 2 higher twist effects become 
important and are not well known [21] . Secondly, the pQCD expansion loses its meaning as 
Q approaches A.qcd where, as a consequence of renormalization and regularization, a™ 
itself is singular. It is necessary to use an appropriate theoretical framework to circumvent 
these difficulties. Such frameworks have been developed, see for example U] and [l^ . We 
use the method of "effective charges" of Grunberg 13] , where the non-perturbative terms and 
higher order perturbative processes are absorbed in the definition of the coupling constant. 
In our case, it obeys the following definition: 

rr n = \g A [l - (3) 
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Eq. 3 provides a definition of an effective QCD running coupling that we will explore here. 
The inherent systematic uncertainties in this experimental Ansatz, and in those of the 
various theoretical approaches are unknown. Their comparison provides a framework for 
further analysis. 

The coupling constant defined with Eq. 3 still obeys the renormalization group equation 
da s (k) /dln(k) = (3(a s (k)) 13]. The first two terms in the Eq. 2 and 3 series are independent 
of the choice of gauge and renormalization scheme. Consequently, the effective coupling 
constant is renormalization scheme and gauge independent, but becomes process-dependent. 
These process-dependent coupling constants can be related by using "commensurate scale 
relations" which connect observables without renormalization scheme or scale ambiguity 1 141 • 
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In this topic, an important relation is the Crewther relation [14J, 1151 ] . 

Considering an effective coupling constant yields other advantages beside renormalization 
scheme/gauge independence: such a procedure improves perturbative expansions JjL 
the effective charge is analytic when crossing quark mass thresholds, is non-singular at 



Q = Aqcd, and is well denned at low Q 2 [16j. The choice of defining the effective charge 
with Eq. 3 has unique advantages: Low Q 2 data exist and near real photon data will be 



available soon from JLab 
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181 ]. Furthermore, sum rules constrain a s>gi at both low and 



p—n 



high Q limits, as will be discussed in the next paragraph. Another advantage is that T± 
is a quantity well suited to be calculated at any Q 2 because of various cancellations that 
simplify calculations 0, l^ . 

Using Eq. [3] and the JLab data, a s>gi can be formed. The elastic (x = 1) contribution is 
excluded in Y\~ n . The resulting a s ^ gi /7i is shown in Fig. [TJ Systematic effects dominate the 
uncertainties, see ref. js] for details. The uncertainty from qa is small. We also used the world 
data of the Bjorken sum evaluated at < Q 2 >=5 GeV 2 to compute a Sjgi /ir. We can also use 



a model for Ti and, using Eq. [3], form a S)Sl . We chose the Burkert-Ioffe model |20j because 
of its good match with the experimental data on moments of spin structure functions [l|- 
. It is interesting to note the behavior of r o 4i9l near Q 2 = where it is constrained by 



the Gerasimov-Drell-Hearn (GDH) sum rule 21[ that predicts the derivative of the Bjorken 
integral with respect to Q 2 at Q 2 

Q 2 



o. 



22|: 



p—n 



lQn 2 a 



(GDH P - GDH T 



-Q 2 



■ K p 



K n ■ 
M 2 ' 



(4) 



where n r 



is the proton (neutron) anomalous magnetic moment and a is the QED 



coupling constant. This, in combination with Eq. 3, yields: 
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The constraint is shown by the dashed line. At Q 2 = the Bjorken sum is zero and a St 
a particular property of the definition of a s>gi . At larger Q 2 , where higher twist effects are 
negligible, a s , 9l can be evaluated by estimating the right hand side of Eq. 2 (using a v ® CD as 
predicted by pQCD) and equating it to g^fl — a SjS1 /7r]/6. The resulting a s>gi is shown by 
the gray band |23j. The width of the band is due to the uncertainty in Aq C d- Finally a s ,F 3 
can be computed using data on the Gross-Llewellyn Smith sum rule [2J], which relates the 
number of valence quarks in the hadron, n v , to the structure function F s (Q 2 ,x) measured 
in neutrino scattering. At leading twist, the GLS sum rule reads: 

2 / . /o2>\ 3 



jf 1 F 3 (Q 2 , x)dx = n v [l- - 3.58 (^^) ~ 20.21 f 



a s (Q 2 
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Using the data taken by the CCFR collaboration [25] , we can apply the same procedure as 
for the Bjorken sum rule to compute a s ^ 3 . We expect a s ^ 3 = a S)gl at large Q 2 , since the 
leading twist Q 2 -dependence of Eq. 2 and Eq. 6 is identical up to , up to a very small 
difference at order coming from the light-by-light contribution to the GLS sum. 

In principle, the commensurate scale relations derived by Brodsky and Lu in ref. [ijj] 
should be applied when comparing a Sigl to a s ,F 3 (and a s>T in Fig. 1). In practice, however, 
the resulting corrections are small; the same scale is used in the a Sjgl and ct s ,F 3 comparison 
and the correction from the light by light contribution decreases the value of a St p 3 by at 
most 1%. These corrections were neglected. The commensurate scale relations should be 
used when comparing a SjT (Q) extracted from OPAL data on r-decay [lij to a s , gl . At leading 
order (the only order available since only one Q value of a SjT (Q) is available) the Q correction 
necessary to compare a S)gl to a s>T leads to a 20% change on the effective charge. Such Q 
correction was applied. 

In summary, assuming the validity of the GDH and Bjorken sum rules for low and large 
Q 2 respectively and using the JLab data at intermediate Q 2 , we can evaluate at SjS1 at any 
value of Q 2 . The absence of divergence in a Sjgl is obvious since it is defined from finite 
experimental data. It is interesting to notice that a s>gl loses its Q 2 -dependence at low Q 2 . 
This feature was suggested by the work of Brodsky et al. based on r-decay data [16] and 
by other theoretical works, as will be discussed below. 
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Many theoretical or phenomenological studies of a s at low-Q 2 are available. See [2a. 127] 
for reviews. The theoretical studies comprise Schwinger-Dyson Equations (SDE), lattice 
QCD, non perturbative QCD vacuum, and analyticity arguments. Phenomenological studies 
are based on quark model spectroscopy, low-x-low-Q 2 reactions, parametrization of nucleon 
and pion form factors, heavy quarkonia decay, and ratio of hadron production cross section 
to production cross section in e + e~ annihilation. 

Similarly to experimental effective charges, different definitions of the strong coupling 
constant at low Q 2 [2^ are possible in the theoretical calculations. How they are related 
is not fully known. Furthermore, these calculations should be viewed as indications of the 
behavior of a s rather than strict predictions. Although some theoretical uncertainties due 
to parameterizations are shown, the existence of unknown systematic/model uncertainties 
should be borne in mind. However, it is interesting to compare the various calculations to 
our result to see whether they show common features. 
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FIG. 1: (color online) a s {Q)/ir obtained from JLab data (up triangles), the GLS sum result from 
the CCFR collaboration [25[] (stars), the world T^~ n data (open square), the Bjorken sum rule 
(gray band) and the Burkert-Ioffe Model. a s (Q)/ir derived using leading order commensurate 
scale relations and the a SiT {Q)/ii from OPAL data is given by the reversed triangle. The dashed 
line is the GDH constraint on the derivative of a Sj91 /7r at Q 2 =0. 

In theory, solving the SDE equations can yield a s . In practice, approximations are 
necessary and choices of approximation lead to different values of a s . In the top panels and 
the bottom left panel of Fig. [2], we compare our data to different approaches 
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331 ] . The uncertainty in Cornwall's result is due to the uncertainties in parameters that 
enter the calculation. The uncertainty on the Bloch curve is due to Aqcd- There is a good 
agreement between the absolute value obtained from the present Ansatz and the results of 



Bloch 



3l| and Fisher et al. |30j], while the results from Maris- Tand y | 32l| and Bhagwat et 



al. 33| do not agree as well. The older calculation from Cornwall [29] disagree with the 
result of our Ansatz. The Godfrey and Isgur curve in the top right panel of Fig. [2] represents 
the coupling constant used in a quark model 34]. 

It is interesting to notice that the Q 2 -dependence of a s>gi and the ones of the calculations 
are similar. A relative comparison reveals that the Q 2 -dependence of the Godfrey-Isgur, 
Cornwall and Fisher et al. results agree well with the data while the curves from Maris- 
Tandy, Bloch et al. and Bhagwat et al. are slightly below the data (by typically one sigma) 
for Q 2 > 0.6 GeV. 

One Ansatz that has received recent theoretical attention for the QCD coupling at low Q 2 
is related to the product of the gluon propagator dressing function with the ghost propagator 
drying taction squared Q. Gluon and ghost propagators have been conrputed in lattice 
QCD by many groups. The lattice results offer a fairly consistent picture [27j and agree 
y with the various SDE propagator results [35(]. It is calculated in particular in 



reasona 
refs. 27 



36] . In the bottom right panel of Fig. El we compare a Si91 to the lattice result from 



ref. [36j. We plot the lattice result that is believed to be closer to the continuum limit 
There is a good agreement between the lattice QCD calculations and our data. 

To conclude, we have formed an effective strong coupling constant a Sj91 at low Q 2 . Data, 
together with sum rules, allow to obtain oi s>gi at any Q 2 . The connection between the 
Bjorken and the GDH sum rules yields a value of a Sj9l equal to n at Q 2 = 0. An important 
feature of a s>gi is its loss of Q 2 -dependence at low Q 2 . We compared our result to other 
coupling constants from different reactions. They agree with each other, although they 
were defined from different processes. We also compared a s>gi to SDE calculations, lattice 
QCD calculations and a coupling constant used in a quark model. Although the relation 
between the various calculations is not well understood, the data and calculations agree in 
most cases especially when only considering Q 2 -dependences. It will be interesting in the 
future to pursue the same analysis with lower Q 2 data that will be available both for the 



neutron [l7j and proton ISfj 
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FIG. 2: (color online) a s ,gi from JLab data and sum rules compared to various calculations: top 
left panel: SDE calculations from Fisher et al. and Cornwall; top right panel: Bloch et al. (SDE) 
and Godfrey-Isgur (quark model); bottom left: Maris-Tandy (SDE) and Bhagwat et al. (SDE); 
bottom right: Furui and Nakajima (lattice QCD). 
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